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Derivation of a General Perturbation Solution: Its Application 
to Determination of Orbit 


An analytical solution to three-body problems has 
been applied to the problem of predicting the orbit 
of a lunar satellite and determining the orbit of a 
near-Earth satellite. Using this solution, the state 
vector may be generated at any desired time without 
intermediate numerical extrapolation. 

General perturbation theory is concerned with 
the development of analytical solutions to the equa- 
tions of motion of a satellite moving under the in- 
fluence of an arbitrary gravitational force. A first- 
order general perturbation solution has been devel- 
oped for the motion of a satellite under the influ- 
ence of an arbitrarily shaped primary body and a 
point-mass third body. A set of nonsingular orbit 
elements is used to describe the satellite’s motion; 
thus the results are valid for all circular and 
elliptic motion. The solutions utilize the development 
of the disturbing function in terms of the Keplerian 
elements as given by Kaula. They are obtained by an 
extension of Ingram’s work to include the third- 
body effects. 

The general perturbation solutions are then incor- 
porated into an orbit-determination scheme. The re- 
sult is a rapid, accurate program that is valuable 
for preliminary design work. Its value as a research 
tool is demonstrated by its use in studying several 
problems associated with orbit determination. 

A modified first-order general perturbation solu- 
tion is developed for the motion of a point-mass 
satellite under the influence of an arbitrarily shaped 
primary body and a point-mass third body. Results 
yielded by these solutions are compared with numeri- 
cally integrated trajectories for various orbital con- 


figurations. It is shown that the analytical solu- 
tions are accurate and are amenable to computer 
evaluation. 

Linear estimation theory is reviewed, and a tech- 
nique for recursive estimation of the observation- 
error covariance matrix is developed. The analytical 
solutions are used as the basis for an orbit-deter- 
mination scheme that compares favorably with an 
orbit-determination program using numerical integra- 
tion. The accuracy of transition matrices, formed by 
numeric partial differentiation, is studied briefly for 
justification of their use in the analytical program. 
The effects of the coordinates used in the estimation 
process are compared numerically, using both 
Keplerian elements and Cartesian coordinates to deter- 
mine the orbit of a near-Earth satellite. In addition, 
the recursive scheme for estimating the radar 
covariance matrix is tested using the same sample 
problem. 

For the numerical studies, the problem of esti- 
mating the orbit of a, near-Earth satellite under the 
influence of the gravitational harmonics J ! 2o> «^30> 
and was considered. Data on range, azimuth, 
and elevation, as observed from four tracking sta- 
tions, were simulated and corrupted with normally 
distributed random noise. 
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